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ABSTRACT 


Stationary cyclotron waves attached to the bow shock are produced 
by the current in the shock. Waves upstream from the shock are circularly 
polarized electron waves, the polarization being in the sense of a 
right-hand screw when the interplanetary field is directed toward the 
shock, and left-hand when directed away. The downstream wave field has 
two components, one from the electron cyclotron branch and one from the 
ion cyclotron branch of the dispersion relation. The latter has a longer 
wavelength and a larger amplitude (by the ion/electron mass ratio) than 
the electron wave. It is elliptically polarized rather than circularly, 
and always in the opposite sense to the upstream electron wave. Therefore, 
the sense of polarization should be seen to reverse as the shock is traversed 
in either direction. When the ambient magnetic field is not perpendicular 
to the shock, there also exists a stationary electrostatic wave with its 
electric field normal to the shock. Both up -and downstream the ratio 
of the electrostatic wave to the magnetic wave amplitudes is of the order 
u/c tan 0, where 0 is the angle between the normal to the shock and the 
magnetic field, and u is the plasma flow velocity. 
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I. Introduction 

x A thin current sheet was used in a previous paper [Tidman and Northrop, 

1968] as a model for the bow shock and as a source of waves extending 
upstream and downstream from it. Ar. example of these waves is reproduced 
in Fig. (l) from Heppner et al. [1966]. So long as the wavelengths are 
much larger than the thickness of the collisionless shock transition, 
it is unnecessary to examine the detailed mechanism of the shock that 
results in the sheet current. In Fig. (l) the bow shock is crossed at 
2 min. 10 sec., where the amplitude and wavelength are seen to change 
suddenly. 

Observed wavelengths [Fredricks, I96S] do not in fact always satisfy 
this condition, and work reported here must then be extended to be applicable. 
This present work makes predictions about the numerous shock crossings 
where the model does apply. 

We use the same calculation method as in Tidman and Northrop [1968] 
tut carry it further, in that we perform the integrations and obtain 
explicit expressions for the wave electric and magnetic fields. There 
only wavelengths were estimated and roughly compared to observations. 

We can now study the amplitude of these cyclotron waves, their 
ellipticity, if any, and the polari7ation sense for comparison vith 
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observations. In addition we have considered a more general geometry 
than there by allowing the upstream and downstream flow vectors u x and 
and unperturbed magnetic fields B x and B 2 to be oblique to the shock. 

Thus the model applies everywhere in the ecliptic plane. Out the ecliptic 
plane is not covered because the flows and magnetic fields have been 
taken as coplanar with the normal to the shock. Extension to the most 
general case is not difficult, but we believe that the work reported 
carries this model as far as is reasonably worthwhile until comparison 
with experiment demonstrates its usefulness. 

We consider only cases where the current j in the shock front is steady 
and normal to the ecliptic plane, as in Figure 2, where the x-axis is 
normal to the shock and the y-axis is normal to the ecliptic. 

Tidman and Northrop [1968] gave consideration to the case where j 
is not steady, but we believe that extending the study to such cases vith 
this model is again more appropriately postponed until comparison with 
observation indicates a need. Nevertheless, we give the formulation for 
non-steady currents. 

In the frame of reference of the flowing plasma, the shock appears 

— 4 

as a current sheet being moved in direction — u x ,2 through the plasma. 
Because the magnetic field of the current j interacts with the plasma 
particles, the shock is an impediment to the plasma flow. Waves can be 
expected, just as in the case of a boat moving through the water. We 
can find expressions for the plasma waves by standard methods of analysis 
applied to the Vlasov -Maxwell system of equations, this being the valid 
set of equations for a collisionless plasma. This set is linearized 
about unperturbed values. For example S (upstream) = 3 X 4 B x ^ 


where 
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Bi^ is the wave field. The wave amplitude is often observed to be 
comparable with Bj., so that a linearised treatment may not always be 
quantitative. But the qualitative predictions \e can make are likely 
to be correct. 

We do the problem as an initial value problem in time, the sheet 
current being ”ero for t < 0. We also obtain the asymptotic form of the 
field as t - 

It is easiest to obtain the oblique flow case (Y ^ 0) from the normal 
flow (Y = 0) case via a Lorentz transformation to a frame sliding parallel 
to the shock. We therefore first obtain the wave fields in the shock 
rest frame for Yi, 2 - 0, 0i, 2 ^ 0* lb is easiest in turn to get these 
by first calculating the fields in the frame flowing with the plasma, 
where the plasma is at rest, and then making a Lorent- transformation to 
the shock rest frame . All of our results are presented in the shock frame, 
and are therefore time independent, the current having been assumed steady 
and the t -• ® limit used. A spacecraft moving through this stationary, 
n oniime -dependent wave pattern will of course see time -dependent fields, 
as in Fig. (l) . 

II. Upstream fields for flow normal to shock and magnetic field oblicue 
(ion mass taken as infinite) 

This case is applicable to the subsolar point. By analysis similar 
to the Appendix of Tidman and Northrop [1968], we find the magnetic and 
electric fields in the shock frame to be: 
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The sheet current has been taken as £ j n e n 6(x) so as to pem.it 

^ n--“ 

arbitrary time dependence; K (k, id) is the plasma dielectric tensor; 
3 


3^ is a unit dyadic perpendicular to k; the vertical bars indicate the 
=* 

determinant; and C is the transpose of the cofactor of the matrix 
-* i^c 2 ^ x 

(K - — -r- I.J; x is a unit vector along the x-axis. The elements of K 
or t 

are rather complicated, particularly when there is a magnetic field. 

(K - k ^5 - 1^) is the same as (- ^ ) in Montgomery and Tidman L 1964, p. 1391 

and is also equal to the matrix given in Eq. (1-20) of Stix [1962, p. ll]. 

The id, ky,and k ? integrations have already been performed in inverting 

the transforms, leaving only the k x integrations in (l) and (2). The k 

vector of the waves is normal to the shock. The plasma has been assumed 

k 2 c 2 ^ 1 

stable, so that all zeroes of w 2 |K - — g — I^,| are on or below the real 


axis in the id plane. There will actually be damping due to finite 
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temperature in a stable plasma, so that at large times the damped solutions 
of uu 2 |K - “i” I,p| = 0 contribute nothing to the uu integration. We 
are assuming that the damping time of the least damped pole is much less 


than the smallest 1 /a> n , or m other words, that fluctuations in the 
sheet current occur in a time scale slow compared to the damping time. 

This assumption could need modification if high enough frequency 
fluctuations are present in the current sheet. The only pole that 
contributes to the uu integral at large times is an undamped one due to 
the plasma flow, it appears as a factor (uu + uu n + u x ) in the denominator 
of the uu integrand. 

To make this paper reasonably self-contained, we now reiterate a few 

salient points from Tidman and Northrop [ 1968 ]. The elements of K for 

a cold plasma are algebraic and therefore uu 2 |K(k,uu) - 5plk„«k -0 

n * 

is of the form f (k x ,uu )7T^[uu - uu^(k x )], where f(k x , uu) does not vanish 

for real k x and uu. The uu^(k x )'s are real. Fig. 3 is a sketch of ^(k ). 

To simulate damping in this formally undamped plasma, uu^(k x ) is replaced 

by uu^Cl^) - ie . The ie is essential to determine which poles of the 

integrand in the k x plane contribute to upstream and which to downstream 

waves. Poles of the k x integrand are found by solving - uu n -k x ui-uu^(k x )+ie = 0. 

ic 

Solutions are k = k + — . am. , where k is a solution for e = 0. 
x ox ui+ ^i 7 ox 

dk ox 

The k x integration for x < 0 - that is, for upstream wa/es, is performed 

by closing the contour at infinity in the lower half k x plane, thus 

picking up contributions from any pole in the lower half plane. From 

dujj 

Fig. 3, duuj /dk < 0; then because u.. + ■ must be negative for a 

* OX CliCq2£ 
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solution k to contribute to upstream waves, ui must be < | | to get 
ox ax QX 

upstream waves. This makes sense physically: the group velocity (vhich is 

negative and therefore directed upstream) must be fast enough to beat the 

flow and allow energy to propagate upstream from the source. Similarly, 

the k^ integration contour for downstream waves closes in the upper half k^ 

,du>. . 

plane and picks up contributions from k QX ' s for which u 2 > | ^ * - |. 

The flow velocity downstream must be sufficiently large to sweep energy 
away from the shock in the downstream direction. 

Considerable analysis is needed in connection with both the determinant 
and the C in Eqn's (l) and (2). The elements of the determinant are given 
in (4) and (5) in a coordinate system with B along the z-axis, whereas in 
Fig. 2, B is oblique to the z-axis. Therefore a rotation of the matrix 

=* c 2 zj 

(K - ■ 1^) must be made. [See Appendix for details.] A second 

complication is that when B and k are oblique, longitudinal and transverse 
modes do not separate and factoring the determinant into the form 
w [uu-uL'^k^)] requires the solution of an algebraic equation of high degree. 
Fortunately, if one expands in the ratio Q e /uu ( of gyro to plasma frequencies, 
this problem is eliminated. For the solar wind fl /uu e ~ 10 3 upstream, and 
therefore the approximation should be excellent. 

By use of Stix’s [ 1962 ] equations (1-21) and(2-42) we have 
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where 


A = (1 . Qei r 1 (\ - £e£ _ 


U)2 


qjp 

U)2 


2n 2 


U)4 


cos 2 Qj; 


( 4 ) 


is the electron gyro frequency and uu^ is the electron plasma frequency. 

In Eq. 's (3) and (4) the effect of ion motion on the waves has been 
neglected (m^ = co ), because, as can be seen in Fig. 3, the upstream waves 
arise only from intersections with the electron branch when the flow velocity 
is greater than the Alfven velocity. The finite ion mass does have a 
small effect on the electron branch, but this we ignore. When we consider 
downstream waves, the ion motion must be retained in order to obtain 
wavelengths longer than upstream, as in Fig. 2. From Fig. 3, uj.(k ) > 

OX ° 

2 2 o 2 

so that » 1. Thus in the approximation e /ou e 2 « 1 


A ° ri ' cos2 9 


(5) 


In the radicals in the curly-bracketed factors of Eq. (3), the first term 
sin 4 0 is dropped as small. This approximation is good so long as 9 is 
not too close to 11/2, at which angle the second term is identically zero. 
However, for the solar wind "too close" is within about 10 radians of II/2, 
so that the approximation is excellent in any practical case. To order 
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, Eq. (3) becomes 




(Q_2_aj2 )a)6 i_ 


k 2 c 2 (uu-fi cos0 )+uxu e 2 | a j k 2 c 2 (urK2 e cos0 )+uxu 2 j. 


where in addition 1 has been dropped compared to ■ , which equals —g . 

For a steady sheet current, j = 0 for n / 0, and uu 0 = 0, so that 

j = j only. The quantity x x C • j Q then appearing in the numerator of 

Eq. (l) is given in the Appendix. In the ^ « 1 approximation it becomes 

^e 

(with j Q = y J y ) 


“ x C • *1 ~ H!fi£ | r ; ^c 2 ^ 2 -^ 2 cos 2 8 ) + nipV + “ iujp 2 Qe cos 0 . 

x J 0 IU2 L U)2 (0)2 - P-2J J U) (o)2 - 0 2 J J J y 


where again 1 has been dropped compared to k 2 c 2/ u) 2 . Similarly in Eq. (2) 
(see Appendix) 


k u.xic • C • = 

v 1 O 


~ ikyUiUUe 2 Qp (u>q 2 + k 2 c 2 )sin 0 
X (1)3 («B2 - n*) 


The ?eros of the right hand side of (6) are needed, with u) set equal 
to - k x u x + ie, rather than the - kj^ indicated in Equations (l) and (2), 
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n 


1^,1 is of the form f(k 
where uu^(k^) must be replaced by uu^k^) - ie to include the damping and 
a) is to be replaced by - ^ x u i* Thus the damping can be equivalently introduced 
by writing uu as - k Uj. + ie and not modifying the uu. (k ). It is then 

X 1 X 

unnecessary to write the determinant in explicitly factored form. 

We thexefore set uu = - k x u i + ie in (6) and solve for k^. (The subscript 
x on k^ will be dropped from here on. ) The first square bracket ir 
Eq. (6) has zeros at k* and k 2 , the second at k 3 and k 4 , where 


for the following reason: K - 


2„2 


kc 

UU2 



fiUcos 9 L , /„ / Ti. 2 cos 2 0 

l = - + 1/2 - 


1/2 


2 u a 


Ul" 


=F ie 


1 /p 

k __ cos 9 _ -j_/ 2 ,' U-. 2 cos 2 9 4uy 2 "• 


2u, 


ur= 


± ie 


( 9 ) 


T /p 

. £e£os_e + v / t<g 2 cos g 8 . Jaa^ y ± lt 
3 2u x U X 2 C 2 ./ 


1^2 


cos 9 , /_ /Q g cos 2 9 _ >, =F ie 

4 “ 2Ux " Uj. 2 c2 y 


The upper signs on the ie terms apply when - H < 9 < and the lower for 
II < 9 < . Very small regions around U/2 and 5n/2 are excluded by our 

approximations, as stated previously. Some positive definite factors multiplying 
the ie terms have been omitted because they do not affect the overall sign 
of the term. We see that k x and k 4 give poles of the integrands in (l) and (2) 
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contributing to upstream waves, while k 2 6111(1 k 3 contribute to downstream 
(if Ui is replaced by u 2 ) for - ^ < 0 < -^ . In other words, uj_ + 

o 

is negative at k x and k 4 . 

The roots k x and k 2 merge and become complex when ^ ~'~ 0 2 ' ~ " S oes 

negative. In Fig. 3 this happens when the line uu = - ku x through the 
origin becomes tangent to the electron branch. This occurs if the solar 
wind velocity u^. becomes too great, greater than about TOO km^sec for 0=0 
(using the "standard" values for the solar wind given in the Appendix). It 
can also occur if the Ui remains fixed and the angle 0 between Bx and the shock 
normal increases beyond a critical value (about 55° by our standard values). 
Beyond the critical values the radical is imaginary and the waves are heavily 
damped, in just a few kilometers, in fact. Our standard upstream values 
make r- = Xi = 52 km at 0 =0. 

It should be noted that in the f e yu> e « 1 expansion, the dispersion 

relation for the transverse waves (square brackets in (6)) for 0 / 0 

are obtained from those for 0 = 0 by changing Qg to Qg cos 0. One can 

imagine in Fig. 3, as 0 is increased from zero, that the electron branch is 

reduced by cos 0, so that eventually it becomes tangent to ui = - ku x . (it 

k 2 c 2 =* \ 

is not, however, true that the elements of the matrix V K - I T ! for 

0^0 can be obtained from the 0=0 case by the same prescription. ) 

Both factors of A in (4) are independent of k and therefore represent 
non-dispersive waves, with zero group velocity, which cannot contribute 


to upstream waves. 
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We give no details of the k^ integrations, which are straightforward and 
require evaluation of residues at k x and k 4 for - — < 0 < ^ , and at k 2 and k 3 
for ~ < 0 < ^ . The explicit expressions (9) for k x and k 2 must be used 
after performing the integrations in order to show that the coefficients 
obtained for the y and z magnetic field components are in fact equal, 
although they appear quite different. Their equality establishes that 
the magnetic field is circularly polarized for all values of 0. 

We find the following for the upstream fields: 

_ n Q n 

For - g < 6 < 2 , 



(10 



4n Ul 2 (uu e 2 + k! 2 c 2 ) J v tan 6 x sin k x x 
k 1 c 3 (Q e 2 c 2 cos 2 9 - lwu e 2 Ui 2 )V 2 


(11 


and for ^ < 0 < ^ , 

^(l), x 4IIu 1 iUp 2 I 1v(y sin k 2 x + z cos k 2 x) 
B (x) = Q - 


(12 


-♦(1), v 4nui 2 (o) e 2 + k 2 2 c 2 )j v tan 0 x sin k 2 x 
E (x) " (-k^ J (ne 2 c2 ^trn^jTT7s— 


(13 


Eq. (9) gives k x and k 2 . The ie can now be omitted 
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Several important conclusions are to be drawn from these expressions: 

1. An electrostatic field is present when the magnetic field is oblique to 
the shock normal (0 / 0). The ratio of the electric to the magnetic 
field amplitude is Ui(uu e 2 + ki 2 c 2 ) tan 9/u) e 2 c. Because increases 
with 9, the maximum value of the ratio occurs at the maximum 9 — i.e.,9 . 
Our standard values make this maximum ratio about 3 Ux/c. Thus the 
ratio will always be < 3ui/c upstream. 

2. The magnetic field is circularly polarized in a right-hand screw 
sense about the component of unperturbed field normal to the shock, 
when that component is directed towards the shock. It is left-hand 
polarized when the field is directed generally away from the shock. 

The B^ is illustrated in Fig. 4 for both cases. 

3. Reversing the unperturbed magnetic field (i.e. changing 0 to II + 0) 
reverses the sense of the polarization, as just stated, but the 
amplitude remains unchanged. 

4. The fields become infinite at either the critical angle or critical 
flow velocity - that is, when (fi e 2 c 2 cos 2 0 - lw) e 2 ui 2 ) = 0* The 
physical interpretation is that at the critical condition the group 
velocity equals the flow velocity so that wave energy radiated by the 
bow shock cannot move away from it and therefore builds up in amplitude. 
Of course, the result is not to be taken too seriously, because 
linearization will become invalid. However, one might expect to see 
larger fields when conditions are near critical. 
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III. Upstream fields for flow and magnetic field both oblique (ion 
mass taken infinite ) 

At a point in the ecliptic plane not at the subsolar pointy the plasma 
flow and magnetic field are both oblique to the shock, as in Fig. 2 . If 
there is a current sheet only in the shock frame, there will in general be 
both a current sheet and a change sheet seen in other frames, in particular 
in the frame in which the plasma is at rest. This charge is of order ui/c 
and is obtained from the charge -current four vector [Jackteon, 1962, p. 578 ]. 

The entire problem could be redone, including now the test charge sheet, which 
one expects will excite electrostatic waves. Fortunately, the order Uj/c 
charge eventually has only an order Ui 2 /c 2 effect on the fields. Because we 
are not working to this order, the charge can be ignored. Moreover, for 
the case we consider, where j is normal to the ecliptic and therefore 
perpendicular to ui, there is no charge to order u x /c even. 

The easiest way to obtain the fields for the oblique flow case is to make 
a Lorentz transformation in the ecliptic plane in a direction parallel to the 
shock plane; the oblique flow then becomes a normal flow, for which we already 
have the fields. Because there is no charge in the shock frame, the current 
in any other frame is the same through order Ux/c, so that j can be taken as 
unchanged when making the Lorentz transformation. The situation in the two 
frames of reference is as in Fig. 5, where the plane of the page is the ecliptic 
plane. The magnitude and angle of the magnetic field to the normal have been 
indicated by the primes as being different in the two frames. There are 
unperturbed electric fields In Figs. 2 and 5 normal to the ecliptic that 
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permit the plasma to flow at an angle to the magnetic field. When there is 
an electric field in one frame, the magnetic field in another frame differs 
in order Ui/c from that in the original frame [Jackson, 1962, p. 300]. 

However, the electric field Ex in the shock frame is (ui/c) Bx sin (0! - Yi) 
into the page. Because it is already of order u^/c, its effect on B x 
is of order u^/c 2 and negligible. So ©i' = 9i and Bj.' = Bj.. (The 
electric field Ei' = (ui/c)Bi sin 0 cos (3. We will not need to use it, 
however. ) 

Equations (10)-(13) give us the fields in the normal flow frame of 
Fig. 5, provided we replace Ui in those equations by the applicable normal 
flow velocity here - namely, Uj_ cos Yj.. We then make a Lorentz transformation 
to a frame moving at u* sin Yi in the negative z direction, as shown, and the 
resulting fields are those seen from the shock frame of Fig. 5 . The 
transformation is straightforward, the vector form given for example in 
Jackson [ 1962 , p. 380] being easy to use. The relations between the wave 
fields in the two frames of Fig. 5 are found to be: 


E 


(1) 


- r_ ( 1 )' . “i sln „ U) r 

= X I " 1 -D 


(14) 



;?(l ) 1 . A ui sin Yi 
= B w f y 1 



( 15 ) 


where the prime means in the normal flow frame and the subscript specifies 
the component. From equations (10) and (ll) (for - II/2 < 0 < II/ 2 ) 


T W 
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(l)' - WlJyUx 2 cos 2 Yi(iu e 2 + ki 2 c 2 )tan 0i sin kix 

" x ki’cs (n e 2c2 cos2 0 i “ H? 2 u i 2 cos 2 Yi)i/2 


( 16 ) 



WlJyUjcos YiOJ e 2 (y sin k i x + z cos k i x ) 

k£ 2 (fi e 2 c2 C0S2 0! - lWJ e 2 Ui2 COS^YiW 2 


(17) 


and similarly from (12) and (13) for n/2 < 0 < 311 ^2 . Note that 
-(1)' , -(1)' 

E is of order n 1 /c smaller than B , so that the second term on the 
right hand side of (15 ) may be dropped and = B^ . The wave fields 

for the oblique flow case are then: 

For - n/2 < 0 < n/2 



ldljyuicos YiU) e 2 (y sin kix + z cos kj,x) 
k^ 2 (Q^c* cos^0 1 - lwu e 2 Ui 2 cos 2 Yi) 1 ^ 2 


(18) 


:(D 


■4IIJ Ul 2 cos Yi[u) e 2 sin (0 x -Yi) + k x 2 c 2 cos Yi sin Oijx sin k x x 
kx c 3 cos 0i (n 2 c 2 cos 2 0! - lwo 2 Ui 2 cosSyj^) 1 /^ 


(19) 


and for n/2 < 0 < 311/2 


—4 

B 


( 1 ) 


4n j v Ul cos Yi u) e 2 (y sin k s» x + * cos 

k 2 C 2 (n e 2 c2 cos 2 0! - 4uJ e 2 Ul 2 cos 2 Yi) 1 / 2 


( 20 ) 
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-( 1 ) 4nj v u x 2 cos Y i [g) e 2 sin (8^-Yi) + k 2 2 c 2 cos Yi s in 0i]x sin k 2 x 
(- k^: 3 cos 0! (Q e 2 c 2 cos 2 9j. - 4a> e 2 U! 2 cos 2 Yi) 1 / 2 


( 21 ) 


where k 1>2 = - 


flp cos 0 1 


f(l e 2 cos 2 9i wa \ 


: + l/p - o ~ 1 

2ux cos Yi ' ' Ui^ cos^ Yx c 2 


1/2 


, the upper sign 


giving kx and the lower ki. 

The magnetic fields (18) and (20) for oblique flow can be obtained from 
the normal flow case (10) and (12) by changing u x to u* cos Y. There is no 
corresponding simple prescription for the electric fields. The polarization 
of the wave remains circular with oblique flow. Indeed if Yi = 0i or 
II + 9 so that the flow and unperturbed field vectors are parallel or 
antiparallel, the magnetic wave fields are independent of Y and 0; for if 
cos Yi = cos 9x, (18) becomes independent of the angles, as does k]_. Similarly 
for (20) if cos Yi = - cos 

The critical condition for oblique flow is reached when the radical in 
the denominators vanishes. As one moves along the bow shock from the subsolar 
point towards the dawn or dusk meridian, the angle between the solar wind 
and the magnetic field remains fixed (at about 4-5°) while their angles to the 
normal changes. For the standard solar wind parameters, we find that there 
should be upstream waves from the subsolar point to the dawn meridian, but 
only for about 10° away from the subsolar point towards th. dusk meridian. 

The magnetic field is not steady at 45° to the solar wind, but fluctuates in 
magnitude and direction. Therefore the 10° limit is not to be regarded too 
rigorously - only the qualitative prediction that upstream stationary waves 
should be less frequently observed on the dusk side. 
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IV. Downstream fields for flow normal to shock and magnetic field oblique 

The downstream case is more difficult than the upstream to treat 
realistically because the plasma is no longer cold. The kinetic energy of 
flow has been converted to ther.al energy by the shock. The conservation of 
mass, momentum, and energy equations relate the downstream conditions to the 
upstream. Nevertheless we can analyse the downstream waves independently of 
the upstream (realizing that the downstream unperturbed quantities are related 
to those upstream) by taking the downstream conditions as' given. We also use 
the cold plasma dispersion dielectric tensor downstream with the belief that 
the principal qualitative result is valid. A warm plasma approximation could 
of course be used, but this again we believe ought to wait until need for more 
refined theory should become apparent from experimental observations. 

In Fig. 3, the electron and ion branches are tangent at the origin 
when 0=0, their common slope being the Alfven velocity. In the infinite ion 
mass approximation of Sections II and III, the Alfven velocity is zero, so that 
the ion branch disappears into the horizontal axis, and the slope of the electron 
branch becomes zero at the origin. The downstream intersection with the 
electron branch will clearly give a shorter wavelength than the upstream, because 
|k 2 | > |ki|. This is at variance with Fig. 1. Indeed for our standard 
downstream conditions, the downstream wavelength would be 0.5 km, while the 
upstream wavelength is 52 km. The downstream electron wave is (for 

n n v 

- ^ < 9 <r ”5 and ion mass = ®). 
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A 

^JyU2^ e 2 (y g i n *2* + z cos k 2 x) 

(- k ^: 2 (0g 2 c 2 cos 2 0 2 - 4<J0 e 2 u 2 2 ) V 2 


( 22 ) 


where k 2 is as in (9), but with downstream parameters, such as the downstream 
plasma frequency, etc. This downstream wave may be present, hut clearly 
is not being seen in magnetometer records such as Fig. 1. The only way to 
get a longer wavelength downstream is from the ion branch. In order for the 
line a) = - ku 2 to intersect the ion branch, u 2 must be less than the 
downstream Alfven velocity (when 0 2 = 0), or as we will see, u 2 must be 
< v | cos 0| in general, where v = Q^c/ujj. 

n n 

Expressions (l) and (2) remain valid for downstream fields, if downstream 

quantities are used. We again find that approximating the elements of the 
^ k 2 c 2 

matrix (K - ^ 1^,) is essential to make the algebra tractable. We take 

the limit m e - 0. This limit does not eliminate the electrons from the 

picture; they still carry current. What it does is move the horizontal 

asymptote of the electron branch at to infinity in Fig. 3, so that the 

line w = - ku 2 would never intersect the electron branch. 

The oblique B 2 necessitates rotation of the matrix again, just as in 

Section II, in order to evaluate the numerator of (l). We make the 

approximations that 0 /u> e and %/u)^ are small, being about 10 2 and 10 3 

1 k 2 c 2 ^ 1 

respectively downstream. The determinant |K - — 2 ■ I_| becomes 



N#ui" " • '.■mv: - ••;.■«••. W"' '• i.’ aJ r\ ' 
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i* - *£ "g - 


U)6 


in 2 f" 4 4“1 

) I k 4 c 4 (u/ 2 -ni 2 )co8 2 02+k 2 c 2 u) 2 u) i 2 (l+cos 2 e 2 )- ^2“^ ] ( 23 ) 


If we set (23) equal to 7.ero, with m = - ku 2 + ie, we find k = k D + iek e 


where 


k 2 - + a)l ^ Uz ^ u)j 2 (l + cos 2 9 2 ) 

O Ug ny~C T ~C0p~e 2 C 2 COS2 0 2 


(24) 


c 4 fy 4 cos 2 e 2 - u)^ 4 u 2 4 
e u|f!^ 2 c 4 k Q 2 cos 2 0 2 


k £ is always > 0 when cos 2 0 2 is such that k Q 2 > 0; k Q 2 > 0 when u 2 2 < v A 2 cos 2 0 2 . 
This means that k Q always gives a downstream wave. The square bracket in 
(23) gives the ion and electron branches in the present approximations 
(See Fig. 6). When 0 2 ^ 0, the ion and electron branches are no longer 
tangent . 

Evaluation of the integral is again straightforward, there being two 
poles, at k Q and - k Q . The downstream wave from the ion branch is 
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Eq. ( 25 ) is valid for all 0 2 . The coefficient of cos k Q x can be shown 
to be positive for all 0 2 such that cos 2 0 2 > u 2 2 /v 2 ^ that is, whenever there 
is downstream wave. The polarization about the x-axis is Lhat of a left-hand 
screw when - H/2 < G < II/2 and right-hand otherwise. Therefore the rule of 
thumb for both upstream and downstream waves is: when the magnetic field is 

pointed generally away from the shock, the polarization is left hand about the 
x-axis, and when toward the shock, the polarization is right hand. Because 

—4 

the normal component of B is conserved across the shock, the upstream and 
downstream waves should always be oppositely polarized. 

The wave in ( 25 ) is elliptically polarized. For 9 2 = 0, the ratio of the 
z component to the y component is (3 - u - u|/v 2 A ). The polarization 
is thus very large for u 2 close to and very sensitive to u 2 /'^. The 
downstream waves thus should show rather large ellipticity. 

There is an electrostatic downstream field when 0 2 ^ 0, given by 
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The ratio of this electric field amplitude to the y component (which is 

the smaller component ) of the magnetic field is of the order U S ^ an ^ 2 . 

Finally, the relative amplitude of the downstream electron wave (22) to 

the downstream ion wave can be compared. The electron wave is about /G^ 

smaller than B ^ in (25). 
z 
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Appendix I - Evaluation of C 

=* k 2 c 2 ^ 

The elements of the tensor (K 3^,) are given in Montgomery and 

Tidman [1964] in a coordinate system with the magnetic field along the ?-axis. 
We have the field at an angle (13/2 - 0) to the z-axis. Our tensor C is 
related to the transposed cofactor of Montgomery and Tidman 's tensor, by 
C i j = a mi a n j Mjj, where M^j is Montgomery and Tidman 's transposed cofactor; 
a is the rotation matrix about the y-axis: 


'sin 9 0 cos 0 \ 

0 10 

cos 0 0 sin 9 

The numerators of the integrands in (1) and (2) are easily shown to be: 


/ 


x x C • J n z ^ny c yy + J n z C yz^ y ^ny C zy + ^nz C zz^ 


<V + V 1 U),?, l = i K + V l)(j ny C xy + J nz C xz } 


y % ^ C yy *W + C yz 


+ z iu n (C Z y j n y + C zz Jnz ^ 


Because we take id =0, and because we consider only the case of sheet 
n 

currents normal to the ecliptic, j all vanish and 

nz 
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i X 0 • Jo - J y (i Cyy - y 0 iy ) 

V 1 • C • J 0 = J y V 1 4 C xy 

Therefore only the three elements of C . C , 

xy' yy’ 

We find them to he 


and C are needed, 
zy 


C = M sin 0 + M cos 9 
xy xy ?y 


C = M 

yy yy 


C = - M cos 0 + M sin 0. 
zy xy zy 


Thus in turn we need only M , M t and M , which are easily obtained by 

xy yy zy* J J 

taking the cofactor transposed of - R/uu 2 in Montgomery and Tidraan. 

We find 


M 


- n e 


xy uj(uu2 - 


(1) 2 i-2_2 

(1 - ^-§- sin 2 0) 

( 1)2 ( 1)2 


M 


yy 


= (i . ^ ("I . Jd§i . ^ ^ ^ . Bl ^e 

v (1)2 y V (1)2 0) 2 - (D 2 (l) 2 (o)2 - Qg 2 ) 


M = 1 , U) ^ 2 ^ 2 C 2 

M zy 0 ) (d) 2 - fi 2 ) (D 2 


cos 0 sin 0 
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and finally 


ixc 





k 2 c 2 u ) p 2 Q , 2 sin 2 
uj 2 u) 2 (u) 2 -fle 2 ) (l- 


- i id= 2 0> cos 6 1 * 

- y a){»2 - a>T" J "j 


m v - 2 f . 2 

> 2 fie Sin 9 (1 - tt! ^ 2 ~) 


k ui x x • C • j - 
x 0 


- ik Ui u), 


= - X 


UO (uo 2 - C^ 2 ) 


1 

» e 2 /uo 2 ) J 
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Appendix II - Standard Values 



Upstream 


Downstream 

1 “* 1 


1 1 

1 


|Bi| 

= 3 x 10 5 gauss 

l®2 

1 = 

12 x 10" 5 gauss 

n l 

= 4 cm 3 

n 2 

= 

9 cm 3 

Ui 

= 400 km /sec 

u 2 

= 

100 km/sec 

®e 

= 1.1 x 10 5 sec 1 

w e 

= 

1.7 x 10" sec 1 

COi 

= 2.6 x 10 3 sec 1 

CUi 

= 

3*9 x 10 3 sec" 1 

^e 

= 5-3 x 10 2 sec 1 

n 

e 

= 

2.1 x 10 3 sec" 1 


= 2.9 x 10 _1 sec" 1 

% 

= 

1.1 sec 1 

V A 

= 33 km^sec 

V A 

= 

88 km^ec 
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Figure 1 


Figure 2 
Figure 3 
Figure 4 
Figure 5 
Figure 6 


Figure Captions 

The bow shock with coherent waves of frequencies near 1 cps, 
observed on November 25, 1964. Scales both uncorrected (the 
farthest left) and corrected for the spacecraft field are 
given. The 12-second periodicity is the satellite spin. 
Geometry of the problem 

Dispersion diagram (schematic) for k parallel to B 

Polarization of the upstream wave 

Frames of reference for the oblique flow case 

Downstream dispersion diagram (schematic) in the approximation 

where m -* 
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